Equilibrium Structure and Radial Oscillations of Dark Matter Admixed Neutron Stars 
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In [Leung et al, Phys. Rev. D 84, 107301 (2011)], we presented our results on using a general 
relativistic two-fluid formalism to study the hydrostatic equilibrium configuration of an admixture 
of degenerate dark matter and normal nuclear matter. In this work, we present more analysis to 
complement our previous findings. We study the radial oscillation modes of these compact stars in 
detail. We find that these stars in general have two classes of oscillation modes. For a given total 
mass of the star, the first class of modes is insensitive to the dark-matter particle mass. They also 
reduce properly to the oscillation modes of the corresponding ordinary neutron star, with the same 
total mass, when the mass fraction of dark matter tends to zero. On the other hand, the second 
class of modes is due mainly to the dark-matter fluid. In the small dark-matter mass fraction 
limit, these modes are characterized purely by the oscillations of dark matter, while the normal 
matter is essentially at rest. In the intermediate regime where the mass fractions of the two fluids 
are comparable, the normal matter oscillates with the dark matter due to their coupling through 
gravity. In contrast to the first class of modes, the frequencies of these dark-matter dominated 
modes depend sensitively on the mass of dark-matter particles. 
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I. INTRODUCTION 

Dark matter (DM) has gained more support for its ex- 
istence from observations [l|, such as galactic rotation 
curves , stability of bars in spiral galaxies 0-0] , cos- 
mological structure formation Q , and gravitational lens- 
ing [9(. However, the properties of DM particles, includ- 
ing their mass and interactions, are still largely unknown. 
Recent data from the DAMA and CoGeNT [III ex- 
periments are consistent with detecting light DM parti- 
cles with mass ~ 10 GeV, which are incompatible with 
the null results from CDMS Q2 and XENON On 
the theoretical side, it has been suggested that isospin- 
violating DM may be the key to reconciling these exper- 
imental results 0, [l5[ . 

Despite the uncertainties on DM properties, it is still 
interesting to ask how DM would affect stellar structure 
and whether one could in turn make use of stellar ob- 
jects to put constraints on DM. The role of DM in the 
first generation stars as stellar seeds, together with the 
possibility of DM annihilation as the first phase in stellar 
evolution, are examined in QUlll. The impacts of DM 
on the evolution and formation of main-sequence stars 
are also examined [20l - [2^ |. The accretion and accumu- 
lation of non-self-annihilating DM particles in the cores 
of planets in our solar system and some pulsars are pro- 
posed to cause changes in their orbits |23l - l25| . Similarly, 
the effects of low-mass (<~ 5 GeV) asymmetric DM parti- 
cles on the solar composition, oscillations, and neutrino 
fluxes have been considered recently |26l - t28| . 

Besides main-sequence stars, one may infer DM parti- 
cle properties through their effects on compact stars [29i — 
l3l| . The effects of DM annihilation on the cooling curves 
of compact stars are also studied in LTil . The re- 
sponse of neutron stars (NS) under non-self-annihilating 



DM models, such as asymmetric matter [351 ] and mirror 
matter (3|| have recently been studied. Neutron stars 
with DM cores are inherently two-fluid systems where 
the normal matter (NM) and DM couple essentially only 
through gravity. The technique used in recent studies of 
the structure of these dark-matter admixed neutron stars 
(DANS) is based on an ad hoc separation of the Tolman- 
Oppenheimer-Volkoff (TOV) equation into two different 
sets for the normal and dark components inside the star 
[HI |H, H3] ■ This approach is motivated by the similarity 
of the structure equations between the relativistic and 
Newtonian ones, but it is not derived from first princi- 
ple. In a previous paper [38j], we made use of a gen- 
eral relativistic two-fluid formulation to study the equi- 
librium properties of DANS and proposed the existence 
of a new class of compact stars which are dominated by 
DM. These stars in general have a small NM core with ra- 
dius of a few kilometer embedded in a ten-kilomcter-sized 
DM halo. In the present work, we describe the formu- 
lation of two-fluid DANS in greater detail. We present 
additional analysis to study the equilibrium properties 
and the radial oscillation modes of these stars in detail. 



The outline of this paper is as follows: In Sec. |TT] we 
present the relevant equations used to study the static 
equilibrium structure, moment of inertia and radial os- 
cillation modes of DANS. In Sec. IIIII we study the static 
equilibrium properties of DANS. Sec. IIVI discusses the 
stability of DANS and their radial oscillation modes in 
detail. Finally, we summarize our results in Sec. [V] We 
use units where G = c = 1 unless otherwise noted. 
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II. TWO-FLUID FORMALISM FOR DANS 
A. Static equilibrium models 

The general relativistic two-fluid formalism was devel- 
oped by Carter and his collaborators (see, e.g., [3{|). It 
has been employed in the study of superfluid neutron 
stars (e.g., |40l443| ). where the two fluids are normal 
and superfluid nuclear matter. In this work, we adopt 
the formulation in to study DANS. Here we shall 
briefly summarize the equations we used and refer the 
reader to [!(| for more details. To find the structure of 
a relativistic two-fluid star, one needs Einstein's equa- 
tions G^u = fSitTpv coupled with a matter source. In the 
two-fluid formulation, the matter source is governed by 
the master function A(n 2 , p 2 , x 2 ), which is formed by the 
scalars n 2 — —n a n a , p 2 = —p a p a , and x 2 = —n a p a . 
The four vectors n a and p a are the conserved NM and 
DM number density currents respectively 49] . The mas- 
ter function is a two-fluid analog of the equation of state 
(EOS) and —A is taken to be the thermodynamics energy 
density. 

For a static and spherically symmetric spacetime ds 2 = 
-e^dt 2 + e x ^dr 2 + r 2 (d9 2 + sin 2 6># 2 ), the Einstein 
equations and the fluid equations of motion reduce to 
the following differential equations 
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where the primes refer to derivative with respect to r and 
A 



B = -^—^ C = ~ 2 1^- ( 5 ) 



dx 2 dn 2 ' dp 

The coefficients Aq, Bq, and Cq are given by 
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The generalized pressure ^ in Eq. @ is computed by 



*(n,p) = A + [in + xp, 



(9) 



where fi = Bn + Ap and x = Cp + An. 

Using Eqs. ©-(H, one may calculate numerically the 
structure of a star by choosing suitable boundary condi- 
tions, namely the central densities of the two fluids n c 
and p c , A(0) = A'(0) = and i/(0) = 0. u(0) is fixed 



by matching the solution with the Schwarzschild metric 
at the star surface r = R. The surface of the NM at 
r = i?NM is defined by the condition n(i?NM) = 0. Sim- 
ilarly, we have the condition p(-Rdm) = to define the 
surface of the DM fluid at r — Rum- In practice we 
choose n or p = 10~ a fm~ 3 to define the surface of either 
fluid. It should be noted that in general the two sur- 
faces are different (i.e., i?NM ^ -Rdm)- The star surface 
r = R, where a matching to the Schwarzschild metric is 
performed, is defined to be the larger one of -Rnm and 
-Rdm- The total mass of the star is computed by 



M 



An 



drr 2 A(r), 



(10) 



while the total particle masses (baryonic masses) of NM 
and DM are computed by 

/■-Rnm 

M NM = 47rm n / drr 2 e x/2 n, (11) 
Jo 

M DM = 47rmDM / drr 2 e x/2 p, (12) 
Jo 

where m„ and ttium are the particle masses of NM and 
DM respectively. It should be noted that the sum of the 
baryonic masses Mnm + Mdm is in general different from 
the total mass M, which also includes the contributions 
from gravitational and internal energies. 



B. Choice of EOS 

In the two-fluid formalism, the master function A plays 
the role of the EOS information needed in the structure 
calculation. In this work, we assume that DM couples 
with NM only through gravity. Hence, the master func- 
tion does not depend on the scalar product x 2 = —n a p a 
and is separable in the sense that 



A(n,p) = A NM (n) + A DM (p), 



(13) 



Anm(«) and Adm(p) being the negative of energy den- 
sities of NM and DM at a given number density, re- 
spectively. These assumptions imply that the coefficients 
A = A® = in our study. 

We choose the APR EOS for NM (and BBB2 EOS 
[45| as well for comparison) and ideal degenerate Fermi 
gas EOS for DM. As discussed earlier, DM candidates 
in the mass range of a few GeV are of great interest 
recently. We shall thus consider fermionic DM particles 
in this mass range. 



C. Moment of inertia 

Besides global quantities like gravitational mass and 
radius, it is also interesting to study the moment of iner- 
tia of DANS since it is measurable and plays an impor- 
tant role in the physics of ordinary neutron stars. 
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The moment of inertia of a rotating star in general rel- 
ativity is defined by J = J/Q in the slow rotation limit, 
where J and f2 are respectively the angular momentum 
and angular velocity. Here we follow the formulation de- 
veloped in [4l| to calculate the angular momentum of a 
two- fluid star in the slow rotation limit. First, we need 
to compute the frame-dragging of the star due to its ro- 
tation. In general, DM and NM can rotate with different 
velocities, Q p and fi„ respectively. The frame-dragging 
to is given by [4l| 



1 (r 4 e -(A + ,)/2 L; y _ 167re (A-,)/2 ( ^ o _ Ao)Ln 



where 



: U! — fL, 



(14) 



(15) 



D. Equations for radial oscillations 

To study the stability and radial oscillation modes of 
DANS, we use the set of equations for radial perturba- 
tions of a two- fluid star derived in Sec. IV-B of ^jj. We 
shall discuss our numerical scheme in solving the eigen- 
value problem in detail. 

Assuming DM has no interaction with NM except 
through gravity (in the sense that the coefficient A = 
Aq = ), and adopting the notation S n — rnW n and 
S p = rpWp, the set of equations governing radial oscilla- 
tions is simplified to 



(BS n ) 



,{u-X)/2 



i B S 'nV 



(22) 



Lp — to — • 



(16) 



It should be noted that Eq. (fT4"| is formally identical to 
the equation obtained by Hartle 46] for one-fluid stars 
except for the nonzero source term on the right-hand 
side. In particular, the corotating case (ft n = £l p ) re- 
duces to the one-fluid result. To integrate this equation 
with higher accuracy, one defines a new variable 



T n = (r 



4^-(A+^)/2 



L'n). 



(17) 



We can decompose Eq. (fT4"|) into two first-order equa- 
tions: 



T; i = 16 7 rrV A ^/ 2 (* -Ao)£„+ 



T> - _Lp(A+'^)/2c 
— 4 L °n 

r 4 



(18) 



(19) 



We integrate the two variables T n and L„ from the origin 
to the surface subject to the boundary condition 



L n (R) 



R (dL n 



(20) 



r=R 



With L n and L p calculated above, the total angular 
momentum of the star is given by 



J = — 



8vr 



drT 4 e (A-,)/2 (Monoi?i + XQpQLp y (21) 



When the master function A is independent of the scalar 



product x 



-n a p a (as we assume in this work), the 



two terms in the integral correspond to the NM and DM 
angular momenta (J„ and J p ) respectively [4l[. We can 
thus define the moments of inertia I n = J n /VL n and I p = 
Jp/ttp. I„ and I p depend on the nonrotating background 
quantities of NM and DM separately. The total moment 
of inertia of DANS is then defined by / = /„ + I p . 



o^-")/^(cs p ) = (- 



Ks' p ] 



1 e v l 2 , 



(23) 



where 



H' = 47rre^ 2+x (p 2 C a + n 2 B° - 2* - ^)H 2 



8 ^ e (,+A)/2 



(pC°S'+nB° S' n ), 



87re A / 2 

H 2 = (xSp + HSn)- 



(24) 



(25) 



Note that the quantities W n and W p are related to 
the radial component of the Lagrangian displacement 
of both fluids by <5££ = 



- x / 2 W p e iuJt /r. 



- x l 2 W n e wt jr and 



e "'-w p e 

The radial oscillation modes can be solved by specify- 
ing the boundary conditions at the core, and the cor- 
rect eigenvalues can be obtained by checking whether 
the boundary conditions at the surfaces of both fluids 
are satisfied. In this problem, we have 4 degrees of free- 
dom at the core, namely S n (0), S p (0), S' n (0) and S' p (0). 
Note that the four variables cannot be set arbitrarily, 
otherwise the boundary conditions of both fluids at the 
surfaces cannot be satisfied simultaneously. 

The boundary condition at the surface of each fluid is 
given by the vanishing of the Lagrangian density varia- 
tion of the fluid. The Lagrangian density variations are 
given by 40] 



An 

n 

Ap 

P 



-A/2 



Wn 



WL 



W 



1 

2 #2, 



- -H 2 . 
2 



(26) 



(27) 



Hence, we require the boundary conditions An/n = at 
r = Rnm and Ap/p = at r = Rdm- 

To find the correct boundary conditions at the core, 
we first choose S n to be unity. S' n (0) is chosen accord- 
ing to the Taylor expansion of S n (r) and the regularity 
condition of S n near the center. Namely, we have 



S n (r)=S n (0y + O(r 5 ), 



(28) 



from which we have 5*^(0) = 35'„(0) [similarly, we have 
S' p (0) = 35 p (0)]. S p (0) is chosen such that the boundary 
condition of the inner surface of a DANS model can be 
satisfied. In practice, we first choose a trial eigenvalue uj 
and a trial S p (0) (without loss of generality, we assume 
Rdm < Rnm in this discussion). Then, we integrate up 
to the inner surface of the DANS model. If the boundary 
condition at the inner surface is not satisfied, a new trial 
of S p (0) is used and the integration is repeated. Once a 
correct trial of S p (0) is found to satisfy the boundary con- 
dition at the inner surface, we continue the integration 
up to the outer surface and check whether the boundary 
condition at the outer surface is satisfied. We obtain the 
eigenvalue if the trial w can satisfy the boundary condi- 
tion of the outer surface. 



III. 



STATIC EQUILIBRIUM PROPERTIES OF 
DANS 

A. General Properties of DANS 



In Fig. |1 (a) we show the mass-radius relations of 



DANS for different DM mass fractions defined by 



M 



DM 



M 



NM 



M 



DM 



(29) 



The NM is modeled by the APR EOS and the DM par- 
ticle mass is ttidm = 1 GeV. 

The case e = corresponds to ordinary NS without 
DM. For increasing e, representing DANS with higher 
DM proportion, we observe two results: first, the max- 
imum stable mass decreases and the stars also have 
smaller radii. For example, for e = 0.2, the mass and 
radius of the maximum stable mass configuration are de- 
creased by 35% and 9% respectively, compared to the 
case with e = 0. Second, the M-R curve flattens as e 
increases, such as the one shown for e = 0.8. For inter- 
mediate DM mass fraction, such as e = 0.4, the curves 
are made of two parts: The curve is flat at large R, but at 
small R, it is qualitatively similar to the curve for ordi- 
nary NS. The above patterns suggest that a new class of 
compact stars exists, when e is sufficiently large. We will 
show that these DM dominated stars are qualitatively 
different from ordinary NS in many aspects. 

We plot in Fig. 1(b) the density profiles of the DANS 
models marked as point X in Fig. 1(a) The upper panel 
is for the case e = while the lower one is for the case 
e = 0.8. Although both models have the same total mass 
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FIG. 1: Upper plot (a): Mass-radius relations of DANS for 
different DM mass fraction e. The NM is modeled by the APR 
EOS and the DM particle mass is 1 GeV. Masses are in unit 
of Mq. Lower plot (b): Density profiles of the star models 
at the point X in Fig. 1(a) for an ordinary NS (e = 0, upper 
panel) and a DM dominated star (e = 0.8, lower panel). Both 
star models have mass 0.475 Mq and R = 11.25 km. For 
the DM dominated star (lower panel), Mnm = 0.0981Mq, 
M dm = 0.3923Af Q and R NM = 6.29 km. 



and radius, the mass in the upper one is contributed only 
by NM, while the mass in the lower one is mainly con- 
tributed by DM. For the case e = 0.8, it is seen that a 
small NM core is embedded in a larger DM halo. 

Next we plot in Fig. [5] two distinctive models for the 

EE). 



case e = 0.4 [marked with triangles in Fig. l(a)| . One 
star is chosen from the part of the M-R curve which 
is similar to that of ordinary NS, while another star is 
chosen from the other side. We choose one model (upper 
panel in Fig. [2J| to have M = O.78OM , R = 8.10 km, 
Rdm — 6.60 km and the other (lower panel) to have 
M = 0.332M©, R = 9.42 km, R NM = 8.84 km. For 
the former model (upper panel), we see that the DM is 
engulfed by NM. However, the situation is reversed for 




FIG. 2: Density profiles for two star models with e = 0.4 with 
the APR EOS for NM and ideal degenerate gas EOS for DM 
with DM particle mass = 1 GeV. The solid (dashed) lines are 
for NM (DM). Upper plot: M = 0.780M©, R = 8.10 km, 
Rdm = 6.60 km. Lower plot: M = O.332M , R = 9.42 km, 
Rnm = 8.84 km. 
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FIG. 3: I /MR 2 of DANS against M/R with the APR EOS 
for NM and ideal degenerate Fermi gas EOS for DM, with 
DM particle mass of 1 GeV. The dashed line corresponds to 
the numerical fitting found by Bejger and Haensel [43|. The 
circles correspond to the two stellar models (e = and 0.8) 
at the point X in Fig. 1(a) 



the latter model (lower panel). 

The scaled moment of inertia I = I /MR 2 of DANS 
is plotted as a function of the compactness (M/R, in 
the unit solar mass/km) in Fig. [3J In [13], Bejger and 
Haensel found an approximate universal relation between 
I and compactness (z = (M JM®)/ '(km/ R)). 



I = 



0.1 + 22' 



z < 0.1, 



I =-(1 + 5*), z > 0.1. 



(30) 
(31) 



This formula is shown as the dashed line in Fig. [3l The 
vertical lines (with arrows) at z = 0.05, 0.1 and 0.15 
represent the range of values of / obtained by a large set 
of EOS models which were used to obtain the formula. 
They can be regarded as the error bars of Eqs. (1301) and 
(13"TT) at those values of M/R. The circles in the figure 
correspond to an ordinary neutron star (e = 0) and DM 
dominated star (e = 0.8) at the point X in Fig. 1(a) 
While the scaled moment of inertia of ordinary neutron 
stars can be modeled approximately by Eq. (1301) and 
(13"TT) . Fig. |3] shows that / of DANS depends sensitively on 
the amount of DM. In particular, for the DM dominated 
sequence e = 0.8, the value of I is significantly smaller 
than that allowed for ordinary neutron stars with the 
same compactness. 



But the slope of the linear variation depends on the stiff- 
ness of the NM EOS. The stiffer the NM EOS is, the 
smaller the contraction results. 

In Fig. 4(a)| we plot R/Ro against Mdm, where Rq is 
the radius without DM (i.e., an ordinary NS), for Mnm 
= 1.4, 1.6, and 1.8 M®. We see that the relative changes 
in radii are linear in Mdm- For a fitting in the form 



R 
R< 



= 1 + cliMdm, 



(32) 



we find that m « -0.5 for the APR EOS and <n « -0.7 
for the BBB2 EOS. The larger magnitude of a x for BBB2 
EOS is due to the fact that this EOS is softer than the 
APR EOS. We also see that ai is essentially independent 
of Mnm ■ The effect from DM is no longer linear for large 
AW 

In Fig. |4(b)| we plot I/Iq against Mdm, where Iq is the 
moment of inertia of a NS without DM, for three different 
values of Mnm as above. The change in I is almost the 
same for the three cases: 



I 

— = 1 
h 



a 2 M 



'■l DM, 



(33) 



with a 2 w -0.6 for the APR EOS and a 2 « -0.9 for 
the BBB2 EOS. Again, the value of a 2 does not depend 
strongly on Mnm- 



B. Linear response of DANS 

Here we study the effects of DM in the core on the 
structure of DANS at the small A^dm limit. In the linear 
regime, namely a DANS with small Mdm, the radius and 
moment of inertia of the star vary linearly with AfoM- 



C. Effects of DM particle mass 

As discussed in Sec. U DM candidates in the mass 
range a few GeV is of great interest recently. We shall 
thus study the effects of DM particle mass ttidm in this 
range. We shall compare the equilibrium properties of 
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FIG. 4: Upper plot (a): R/Ro is plotted against Mdm (in M ) 
for the APR and BBB2 EOS for NM and ideal degenerate 
Fermi gas EOS for DM, with DM particle mass of 1 GeV. R 
is the NS radius without DM. Lower plot (b): I/Io against 
Mom (in Mq). Iq is the moment of inertia without DM. 



DANS with different todm- In this part, we use the APR 
EOS to model the NM. 

First, we plot in Fig. 5(a) the total mass M against 
the minimum radius i? m ; n of stable stars allowed for a 
given mDM- We present the results for three different 
cases mDM = 1, 4, and 7 GeV. We also plot the M-R 



curve (solid line) for ordinary NS without DM for com- 
parison. For fixed M and mDM] the radius of the star 
decreases as the mass fraction of DM increases (i.e., the 
star becomes more compact). We define i? m in to be the 
minimum radius below which the star becomes unstable. 



For example, Fig. 5(a) shows that the minimum radius 
allowed for a DANS with total mass M = 1M Q and DM 
particle mass ttidm = 2 GeV is about 10 km. We also 
see that i? m i n increases with m-DM for a given M. 

In Fig. 5(b)| we plot the density profiles of two stars 
with the same M = 1.4M and M DM = O.OLMq, but 
with different DM particle mass tudm = 1 GeV (upper 



FIG. 5: Upper plot (a): Total mass M is plotted against 
the minimum radius R m in of stable stars allowed for a given 
DM particle mass moM- The solid line is the M-R curve for 
ordinary NS without DM for comparison. Lower plot (b): 
Number density profiles of NM (solid lines) and DM (dashed 
lines) for two stars with the same M = \AMq and Mdm = 
O.OIMq, but with 771dm = 1 GeV (upper panel) and 7 GeV 
(lower panel). 



we see 



panel) and 7 GeV (lower panel). In Fig. 5(b) 
that with a higher DM particle mass ytidm = 7 GeV, the 
DM core shrinks to a very small size of about 0.3 km, 
compared to 4.3 km for the case todm = 1 GeV. Also, 
the density of the DM core is much higher than that of 
the NM for m DM = 7 GeV. 

Finally, in Fig. [S] we plot the maximum stable mass 
for DANS models with e = 0.1 but with todm ranging 
from 1 to 5 GeV. The maximum stable mass decreases 
with 771dm- The reason is as follows: the NM and DM are 
assumed to be noninteracting (except through gravity), 
and the DM core is supported only by its own degener- 
ate pressure. It is well known that the maximum mass 
limit for a self-gravitating Fermi gas decreases as the par- 
ticle mass increases. Hence, the onset of the collapse 
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FIG. 6: Maximum stable mass M max is plotted against the 
DM particle mass mx for a fixed amount of DM specified by 

6 = 0.1. 



of a degenerate DM core is responsible for the depen- 
dence of M max on rriDM as seen in Fig. [6j Furthermore, 
the allowed mass fraction of DM inside a stable DANS 
decreases significantly as TOdm increases. For example, 
stable DANS models with Mdm ~ 0.1M© can only be 
formed by DM particles of mass less than about 3 GeV. 



IV. RADIAL OSCILLATIONS OF DANS 

In this section we study the radial oscillation eigen- 
frequencies and eigenfunctions of DANS. We show that 
all DANS with central energy density less than that of 
the maximum mass configuration, regardless of the mass 
fraction e of DM, are stable. We also study the effect of 
e and DM particle mass ttidm on the radial oscillation 
modes. In this section, we use the APR EOS to model 
the NM. Unless otherwise noted, the DM particle mass 
rriDM is chosen to be 1 GeV. 



A. One-fluid limit 

We first present some tests to check the validity of 
the numerical code. We calculate the oscillation modes 
of ordinary NS modeled by the APR EOS in the one- 
fluid limit using our two-fluid code. In practice, this is 
achieved by setting the central density of DM to a suffi- 
ciently small number so that the star is essentially com- 



posed of NM only. In Fig. 7(a) we plot the total mass 



(upper panel) and fundamental mode frequency squared 
(lower panel) against the central energy density. We set 
the central number density of DM to be 8 orders of mag- 
nitude smaller than that of NM in the calculations. As 



expected from the study of ordinary NS, Fig. 7(a) shows 
that the fundamental mode frequency passes through 
zero at the central energy density corresponding to the 
maximum mass configuration. The point u> 2 — marks 
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FIG. 7: Upper plot (a): Total mass (upper panel) and fun- 
damental mode frequency squared (lower panel) are plotted 
against the central energy density for ordinary NS modeled 
by the APR EOS. Lower plot (b): Frequency squared for the 
first three modes are plotted against the central energy den- 
sity for ordinary NS modeled by the APR EOS. The dashed 
lines are obtained by the two-fluid code in the one-fluid limit. 
The solid lines are obtained by a one-fluid code (see text). 



the onset of instability. Beyond this critical density, the 
stars are unstable against radial perturbations. 

We compare in Fig. |7(b)| the frequencies of the first 
three modes calculated separately by the two-fluid code 
(dashed lines) and a different code (solid lines) based on 
the standard one-fluid formulation 481. We see that the 



two sets of mode frequencies agree very well. 



B. Oscillation modes of DANS 



In Fig. 



we plot the total mass (upper panel) and 



the fundamental mode frequency squared (lower panel) 
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FIG. 8: Upper plot (a): Total mass and fundamental mode 
frequency squared are plotted against the central energy den- 
sity for e = 0.2. Lower plot (b): Same as Fig. 8(a) but for 
e = 0.8. 
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FIG. 9: Upper plot (a): The frequency squared for the first 
four oscillation modes are plotted against e. The total mass 
M — 1.4M© is fixed. The circles indicate the absence of 
the corresponding modes at the limit e = 0. Lower plot (b): 
The Lagrangian density variations of the first four modes of 
a DANS with e = 0.005 and total mass M = 1.4M Q . 



against the central energy density for e = 0.2. Similar 
to the ordinary NS (one-fluid) case, the mode frequency 
passes through zero at the central energy density cor- 
responding to the maximum mass configuration. The 
stars are unstable beyond this critical density. For DANS 
with central density lower than the critical density, they 
are stable against radial perturbations. In Fig. |8(b)| we 
show the case e — 0.8 for the DM dominated sequence. 
Our results confirm the stability of DANS. In particular, 
the new class of DM dominated compact stars with a 
NM core embedded in a ten-kilometer sized DM halo are 
shown to be stable. 

In Fig. 9(a) wc plot the first four mode frequency 
squared of DANS as a function of e. The total mass 
M = IAMq is fixed. We can see from the frequency of 
the fundamental mode (n = 1) that increasing e (i.e., the 
mass fraction of DM) has the effect of decreasing the sta- 



bility of the star. The fundamental mode frequency drops 
sharply to zero for e slightly above 0.2, hence indicating 
the onset of instability. We also see that the frequencies 
of the higher order modes (n = 2, 3, 4) in general increase 
with e. It is interesting to notice that the second (n = 2) 
and third {n = 3) modes are missing in the one-fluid limit 
when e = 0. In Fig. 9(b) we show the Lagrangian vari- 
ations of NM [Eq. d2Stjand DM [Eq. ([27])] for the first 
four modes of a star with e = 0.005. The solid (dashed) 
lines are the profiles of NM (DM). Note that the DM core 
extends to about OAR, where R is the radius of the star. 

First let us consider the n — 1 and n — 4 modes. They 
have proper limits at e = 0. For the fundamental mode 
(ri = 1), in the DM core where the two fluids coexist, 
the density variations of NM and DM are in phase. For 
the n = 4 mode, the density variation of NM (DM) is 
larger (smaller) than zero in the DM core. Hence, the 




FIG. 10: The Lagrangian density variations of the first four 
modes for the ordinary DANS with e = 0.1 at the point Y in 
Fig-fTRI 



FIG. 11: The Lagrangian density variations of the first four 
modes for the DM dominated compact star with e = 0.8 at 
the point Y in Fig. |l(a)| 



two fluids are counter-moving in this case. On the other 
hand, for the n = 2 and n = 3 modes, we see that the 
density variation of NM is much smaller than that of DM. 
These modes are dominated by DM motion in the DM 
core. They do not exist when there is no DM (i.e., e = 0). 
However, they emerge even for a very small mass fraction 
of DM. 

To further study the oscillation modes of DANS, we 
choose the two stellar models at the point Y in Fig. 1(a) 
The stars have the same mass M = O.541M0 and radius 
R = 9.68 km, but with different mass fraction of DM. The 
model with e = 0.1 is an ordinary DANS. The model with 
e = 0.8 is a DM dominated compact star with a small NM 
core embedded in a ten-kilometer sized DM halo. The 
Lagrangian variations of the first four oscillation modes 
for the model with e = 0.1 are plotted in Fig. [TO] For 
the n = 1 and n = 4 modes, the general patterns of the 
modes are qualitatively the same as the case e = 0.005 
in Fig. |9(b)| The two fluids are in large part comoving 
(counter-moving) for the n = 1 (n = 4) mode. For the 
n = 2 and n — 3 modes, which are dominated by DM 
motion in the case e = 0.005, we now see that there 
are significant density variations of NM near the stellar 
surface. This can be understood by the fact that the two 
fluids are coupled through gravity. For a very small e, 
the motion of a small amount of DM basically has no 
effect on the NM. Hence, the NM essentially decouples 
from the DM and does not move in the n = 2 and n = 3 
DM dominated modes. However, when the mass fraction 
of DM is comparable to NM, the coupling between the 
two fluids becomes stronger and hence we can see a large 
density variation of NM. 

The Lagrangian variations of the first four oscillation 
modes for the DM dominated model with e = 0.8 are 
plotted in Fig. [IT] We see that for the "ordinary" n = 1 
and n = 4 modes the maximum density variations of the 
two fluids near their surfaces are comparable. Note that 
the surface of the NM core is at about 0.6i?, where the 



radius of the star R is defined by the radius of the DM 
halo in this model. However, for the n = 2 and n = 3 
DM dominated modes, the maximum density variation 
of DM is much larger than that of NM. 



C. Effects of DM particle mass 

To end this section, we study the effects of DM parti- 
cle mass on the mode frequencies. In Fig. fl2l we show 
the frequencies of the first five modes of a DANS with 
M = IAMq and M DM = O.OlAf© as a function of DM 
particle mass Todm- For iiidm = 1 GeV, the second 
and third modes correspond to the n — 2 and n = 3 
DM dominated modes studied above. We see that the 
frequencies of these modes increases with TOdmj while 
the other modes are essentially independent of Todm- 
It is known that the radial oscillation mode frequency 
(squared) scales with the density of the star. Also, for 
a fixed DM core mass, the DM core becomes denser as 
todm increases (see Fig. |5(b)| . Hence, the frequencies of 
the DM dominated modes depend strongly on todm- 

It is also interesting to notice that, while the second 
and third modes are DM dominated modes in the case 
mdm = 1 GeV, this is in general not true for other Todm- 
For example, in the case todm = 3 GeV, it is the second 
and fourth modes that are DM dominated modes. For 
todm = 5 GeV, the DM dominated mode first appears 
only in the third mode. 



V. CONCLUSIONS 

In this paper, we have studied the equilibrium proper- 
ties and radial oscillation modes of DANS using a general 
relativistic two-fluid formalism. We model the NM by 
realistic nuclear matter EOS. The DM particles are as- 
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FIG. 12: Frequency squared for the first five oscillation 
modes are plotted against the DM particle mass ttidm- The 
total mass M = 1.4M© and DM core mass M D m = O.O1M 
of the stars are fixed. 

sumed to be non-self-annihilating and described by an 
ideal degenerate Fermi gas. Our results suggest that 
the structure of these stars depends strongly on the DM 
fluid. In particular, we found a new class of compact 
stars which are dominated by DM. These stars in gen- 
eral have a small NM core with radius a few km embed- 
ded in a larger ten-kilometer-sized DM halo. Since only 
the NM core can emit thermal radiation, the detection 
of a compact star with a thermally radiating surface of 
such a small size could provide a strong evidence for the 
existence of DANS. Furthermore, these DM dominated 
stars also have rather different mass-radius relations and 
scaled moment of inertia compared to ordinary NS with- 
out DM. We have also studied how the radius R and 
moment of inertia I of a star with fixed NM baryonic 
mass Mnm change as the DM core mass Mdm increases. 
In the small Mdm limit, we see that R and / decrease 
linearly as Mdm increases. The slopes of the linear vari- 
ations depend essentially only on the NM EOS, but not 
on the value of Mmm- 

We have performed a radial perturbation analysis and 



studied the oscillation modes of DANS in general. The 
stability of DANS is shown explicitly by calculating the 
frequency of the fundamental mode. For a sequence of 
stars with a fixed DM mass fraction, we see that the 
fundamental mode frequency passes through zero at the 
central energy density corresponding to the maximum 
mass configuration. Similar to the analysis of ordinary 
NS, this point marks the onset of instability. 

Besides the fundamental mode, we have also studied 
the first few higher order oscillation modes. We see that 
DANS in general have two classes of oscillation modes. 
The first class of modes has proper limit when the DM 
mass fraction tends to zero, namely these modes reduce 
to the same set of modes for ordinary NS without DM. On 
the other hand, the second class of modes is due mainly 
to DM. In the limit of a small DM mass fraction, these 
modes are characterized purely by the oscillations of DM. 
The NM fluid is essentially at rest. In the intermediate 
case where the mass fractions of NM and DM are com- 
parable, the NM fluid oscillates with the DM fluid due to 
their coupling through gravity. On the other hand, the 
amplitude of DM oscillations is much larger than that 
of NM in the case of DM dominated stars. We also see 
that the frequencies of these oscillation modes depend 
strongly on the DM particle mass. 

Finally, it should be pointed out that the formation 
mechanism of DANS is not clear. However, our main 
focus in this work is to study the properties of these the- 
oretical objects, if they exist. The formation process of 
these compact dark-matter compact objects cannot yet 
be modeled in current iV-body simulations, which mainly 
focus on the structure formation in galactic or cosmolog- 
ical scales. 
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